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ABSTRACT 

We have developed a new instrument that can meas- 
ure fast transient birefringence and polymer chain ori- 
entation angle in complex fluids. The instrument uses 
a dual-crystal transverse electro-optic modulator with 
the second crystal’s modulation voltage applied 180° 
out of phase from that of the first crystal. In this man- 
ner, the second crystal compensates for the intrinsic 
static birefringence of the first crystal, and it doubles 
the modulation depth. By incorporating a transverse 
electro-optic modulator with two lithium-niobatc 
(LiNb0 3 ) crystals oriented orthogonal to each other 
with a custom-designed optical system, we have pro- 
duced a very small, robust instrument capable of fast 
transient retardation measurements. By measuring the 
sample thickness or optical path length through the 
sample, we can calculate the transient birefringence. 
This system can also measure dichroism. 

We have compared the calibration results and retar- 
dation and orientation angle measurements of this in- 
strument with those of a photoelastic modulator (PEM) 
based system using a quarter wave plate and a high- 
precision 1/16-wave plate to simulate a birefringent 
sample. Transient birefringence measurements on the 
order of 10 9 can be measured using either modulator. 

INTRODUCTION 
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Figure 1. Birefringence Measurement System showing 
polarization and dual-crystal modulator geometry . 

Transverse electro-optic modulators can be made 
very small, rugged and lightweight. LiNb0 3 crystals 
are non -hygroscopic, and they can be mounted in a 
resonant circuit to reduce electrical power needed to 
drive the modulator. These are all highly desirable 
qualities when designing an experiment suitable for 
space flight. 

We have demonstrated through the Mueller matrix 
calculations and by experimentation how to design a 
flight-worthy electro-optic ally modulated birefrin- 
gence measurement system. We show that the retar- 
dation and thus the birefringence of an optically trans- 
parent fluid polymer can be measured accurately using 
this novel technique. 


Existing polarimcters traditionally incorporate cither 
a single-crystal PEM or a rotating retardation plate as 
the modulator clement. Although these systems arc 
capable of very sensitive birefringence and dichroism 
measurements, they are either too large or too fragile 
to survive the forces encountered by space-based ex- 
periments during launch and landing on board the STS 
or sounding rockets. For example, sounding rocket 
launch loads are typically on the order of 1 7g. In order 
to obtain the desired sensitivity while maintaining a 
stable and rugged configuration, we suggest using a 
dual-crystal transverse electro-optic modulator (EOM) 
geometry. The modulator crystals must be oriented 
with their principal axes orthogonal to each other and 
at ±45° with respect to the polarization state of the 
incident light. The drive voltage applied to the second 
crystal must be 1 80° from that of the first crystal. The 
suggested optical configuration is shown in figure 1 . 


OVERVIEW OF BIREFRINGENCE 

Birefringence (An') is the difference in a material’s 
refractive index with respect to the polarization state of 
light propagating through it. In an optically transpar- 
ent stressed polymer material, the orientation and de- 
gree of deformation of the molecular polymer chains 
cause an anisotropic polarization. If the light propa- 
gates through a birefringent material along the z- 
direction, the x and y components of its electric field 
vector will be different from one another. This results 
in a phase difference or retardation (5) as the light 
traverses the sample material. If the material is not 
dichroic, then the retardation is related directly to its 
birefringence. This can be expressed as 


A n„ - n e 
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where n Q and n e are the refractive indices for ordinary 
and extraordinary light rays respectively, X is the vac- 
uum wavelength of the light, and d is the sample mate- 
rial thickness. 

One can apply this knowledge using the stress optic 
law, which provides a simple relationship between the 
material stress and the refractive index tensor. In 
terms of shear stress (x xy ), the stress optic law can be 
expressed as 


T _=— CAn , sin(22’(0) (2) 

2 


where C is the materials stress optic coefficient and 
X(t) is the time-variant orientation angle of the mo- 
lecular polymer chains with respect to the flow or 
stress direction. Thus, in order to measure a material’ s 
shear stress, we must measure the birefringence and 
average molecular chain orientation angle simultane- 
ously. 

THEORY OF DUAL CRYSTAL 
TRANSVERSE ELECTRO-OPTIC 
MODULATION 


We developed a dual-crystal modulation technique 
using the Stokes parameters and Mueller matrix repre- 
sentations of the optical elements. We implemented 
the theory using the transverse electro-optic effect or 
Pockel’s effect and its application to the LiNb 03 crys- 
tals as two separate Mueller matrices. 

A birefringent crystal possesses two distinct indices 
of refraction. These indices are represented by n Q and 
n e , which are the ordinary and extraordinary index of 
refraction, respectively. When a voltage is applied 
across the crystal in a direction normal to the direction 
of light propagation, the difference between polariza- 
tion-sensitive refractive indices changes. This change 
results in a voltage-sensitive phase retardation between 
the ordinary and extraordinary light rays, which can be 
mathematically expressed as 


In . .. 2n 3 VL 

-r (n a -n e )L-— r 33 n 0 — (3) 


where r 33 is the appropriate electro-optic coefficient, L 
is the crystal length, D is the distance between elec- 
trodes and V is the voltage applied across the crystal. 
If V is an alternating voltage, we can modulate the 
electric field vector or polarization state of the light 
passing through the crystal. For pure LiNb0 3 , n 0 = 


2.2866 and n c = 2.2028 at 20°C*. A transverse electro- 
optic crystal configuration is shown in figure 2. 



Figure 2. Transverse electro-optic modulator geome- 
try. Light propagates in the z-direction . We apply 
modulation voltage in the ±y-direction. 

From equation (3), we can see that the electro-optic 
effect is maximized when the electric field vector of 
plane-polarized light is incident on the left crystal face 
at a 45° angle with respect to the z-axis. We can also 
observe that there is static birefringence that depends 
largely on crystal length L. In many crystals, includ- 
ing LiNb0 3 , this static birefringence is highly tem- 
perature dependent 2, 4 , Wc have experimented with 
single-crystal LiNb0 3 modulators and have measured 
the effective static retardation to be 2.24x10 2 radians 
in a 40mm length crystal at 20°C. We see that this 
retardation corresponds to a static birefringence of 
5.69x10 8 radians. Although the crystal’s static bire- 
fringence is quite small, it is not thermally stable". One 
may use a single-crystal transverse electro-optic 
modulator in this type of instrument, but the instru- 
ment must be operated in a temperature-stabilized en- 
vironment. The temperature-dependent static birefrin- 
gence is an undesirable property since it reduces the 
instrument’s sensitivity and stability when measuring 
very small birefringence levels in a sample material. 
Fortunately, by using the dual-crystal configuration 
shown in the preceding figure 1, we can compensate 
for both the residual static birefringence and the ther- 
mal instability while doubling the modulation depth. 

In order to derive how this geometry works, we need 
to describe how polarized light interacts with the opti- 
cal elements. Since light is an electromagnetic wave. 


* From Manufacturer’s data sheet. 
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the Stokes vector, S that consists of four parameters, 
can characterize its polarization state. These parame- 
ters are I, Q, U, and V, which are based on the electric 
field vector. 


r 


\Ex(t)) 2 + (Ey(t)) 2 ' 
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(Ex(t)) 2 -(Ey(t)) 2 
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2Ex(t)Ey(t) cos (S) 
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In equation (4), Ex(t) and Ey(t) are the time- 
dependent x and y electric field vector components, 
and 5 is the instantaneous phase retardation between 
- 7 t and +7t. The parameter I, represents the intensity of 
the light while the parameters Q, U, and V specify its 
state of polarization. 

The polarization state of an optical system can be 
modeled using the Stokes vector of the incident light 
multiplied by the Mueller matrices of the individual 
optical elements. The multiplications are done in re- 
verse order to generate the output Stokes vector that 
yields information about the subject material under 
study. That is, if we represent the incident light with 
the Stokes vector [Sj] and have polarization optical 
elements represented by their 4x4 Mueller matrices 
[Analyzer], [Sample] and [Modulator] as shown in 
figure 3, the output Stokes vector is computed as 

[S 0 J = [AnalyzerhSamplejfModulatorKSd- (5) 

For the dual-crystal polarimeter in our experiment, 
we have the following polarization state configuration 
shown in figure 3. All angles are referenced looking 
into the laser. We define the input polarization at zero 
degrees. The dual-crystal transverse electro-optic 
modulator crystals are at ±45° and the analyzer is at - 
45° with respect to the input polarization state. 
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Figure 3. Polarization state representation of optical 
elements. 


We multiplied the representative Mueller matrices in 
reverse order to calculate an equivalent system matrix 
that was then multiplied by the input Stokes vector. 

The resulting Stokes vector describes the intensity and 
polarization of the light emerging from the last element 
in the system matrix. Once equation (5) is computed 
using the matrix elements, it is Fourier-Bessel ex- 
panded so that I(t) can be written in terms of quantities 
measurable by two heterodyning lock-in amplifiers 
The time-variant intensity, I(t), is given by the first 
element of the resultant Stokes vector and is given 
below. 

/(f) = ldc[ 1 + 2 Jf\) sin(cot) sin(<5 (0) cos(2^ (0) 

+ 2 J 2 (A 0 )cos( 2 tuf){l-cos((y(/))}sin( 2 ^( 0 )cos( 2 /(/)) 
+ ...] 

( 6 ) 

where A 0 = 2.4048. Equation (6) gives us I(t) in terms 
of the in-phase and in-quadrature signal components. 

INSTRUMENT STABILITY 

Since thermal stability is a major problem for a sin- 
gle-crystal LiNb0 3 transverse electro-optic modulator, 
we performed stability tests. We configured the sys- 
tem with the input polarizer at 0°, the single-crystal 
modulator at +45°, and the analyzer at -45°. Wc 
placed a quarter wave plate in the sample space and 
rotated it to give approximately 1/16 wave retardation. 
This is where the I m and I 2(0 signals are equal. We then 
separated these signals slightly by rotating the quarter 
wave plate a few degrees. We increased the laboratory 
temperature by 10°F over a period of 16 hours. Wc 
recorded data over this period to observe the signal 
variance. The stability of the single-crystal modulator 
is very poor due to the temperature dependence of the 
intrinsic birefringence of lithium niobate 4 . We then 
repeated this test using a dual-crystal LiNb0 3 modu- 
lator. Figure 4 shows the results of both tests. 
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Figure 4a. Long term stability test results for a single- 
crystal EOM. The plots show the photodetector and 
lock-in amplifier signal outputs . Idc is the dc signal 
plus the detector bias voltage while I CO and 12 co are 
the first and second harmonic signals respectively. 



Figure 4b. Long term stability test results for a dual- 
crystal EOM . The plots show the photodetector and 
lock-in amplifier signal outputs. Idc is the dc signal 
plus the detector bias voltage while I co andI2co are 
the first and second harmonic signals respectively. 


INSTRUMENT CALIBRATION 


measurement of the dc signal while incrementing the 
modulation voltage for three different retardation val- 
ues in the sample space. The results are shown in fig- 
ure 5a. 



Figure 5a. f, Zero-Bessel” plot to calibrate for 
Jo(A)=0 . In this case, the dc signal is plotted against 
modulation voltage for three different retardation val- 
ues in the sample space. 

We observe that the approximate modulation voltage 
setting is 2V. This calibration method is subject to 
noise in the dc signal and is not as accurate as wc 
would want. 

There are several other ways to calibrate this system, 
but perhaps the most useful method is one that maps 
out the Bessel functions Jj(A) and J 2 (A) simultane- 
ously, yielding a full calibration sweep instead of just a 
single-point calibration measurement. With the instru- 
ment configured so that the input polarizer is at 0°, the 
electro-optic modulator crystals are at ±45°, and the 
analyzer is at -45°, we place a retardation plate in the 
sample space. As long as the retarder is oriented such 
that the first and second harmonic signals are present, 
we may increment the modulation voltage, record the 
in-phase and in-quadrature lock-in amplifier output 
signals, and scale these accordingly to observe their 
correlation with Ji(A) and J 2 (A). The parameter A is 
proportional to the modulation voltage. The calibra- 
tion plot then gives us our modulation voltage at any 
desired Bessel function value. We plotted the correla- 
tion between our electro-optic modulator and the theo- 
retical Bessel function values in figure 5b. 


Since we wanted to set the instrument where J 0 (A)=0, 
we performed a standard "zero-Bessel calibration” as 
described in reference 6. This method involves the 
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Figure 5b. Calibration curves generated using a dual- 
crystal transverse electro-optic modulator. Modula- 
tion voltage was incremented from 0 to 3.5V rms with 
the modulation frequency held constant at 42KHz- 

From the Figure 5b calibration plot, we may more 
accurately observe that Jq(A)= 0 when our modulation 
voltage is approximately 1 .95V rms. The calibration 
plots determine the correct modulation voltage where 
Jo(A)=0. 

By selling the modulation voltage such that J 0 (A)=0, 
we can easily obtain the sample retardation, 5(t) and 
sample orientation angle, %(l). In this case, the value 
of A is 2.4048. The Bessel function calibration values 
are given as Jq(A)=0, J i (A)= 0.5 191 and J2(A)=0.4317. 

MEASUREMENT RESULTS 

We took measurements of the retardation and orien- 
tation angle of a rotating quarter wave plate using the 
eleclro-optically modulated system and the same opti- 
cal system incorporating a photoelastic modulator 
(PEM) in place of the EOM. The intrinsic static strain 
retardation of the dual-crystal EOM was 3.85 x 10 4 
while that of the PEM was an order of magnitude 
greater. We placed a quarter wave retardation plate in 
the sample space. We altered the retardation and ori- 
entation angle using a piezo-driven rotary mount to 
rotate the retardation plate at constant angular velocity. 
The dc signal component as well as the first two har- 
monics were recorded to calculate retardation and ori- 
entation angle simultaneously. These measurement 
results are shown in figure 6. We generated these 
plots using the Extensional Rheology Experiment 
flight prototype optical assembly. The optical assem- 
bly consists of a diode laser head, laser- to- fiber cou- 
pler, polarization-maintaining optical fiber, collimating 


microlens, beam splitter polarizer cube, polar- 
izer/collimator mounting assembly with rotational ad- 
justment capability, beam splitter analyzer cube, ana- 
lyzer housing with rotational adjustment capability, 
analyzer-to-detector mounting fixture, and a high- 
speed silicone pin-diode based photodetector. All of 
these components were used in the optical configura- 
tion to compare performance between the PEM and 
EOM systems. The only component that we changed 
for this comparison was the modulator itself. 


Figure 6b. Retardation (Delta) and Orientation Angle 
(Chi) vs. Time for a rotating quarter wave plate in the 
EOM -based system. 

The figure 6 plots give the retardation ranging from 
zero to nil and the orientation angle ranging from zero 
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Figure 6a. Retardation (Delta) and Orientation Angle 
(Chi) vs. Time for a rotating quarter wave plate in the 
PEM -based system. 
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to rc/8. The slight variation in the orientation angle on 
the figure 6b plot is due to the inability to rotate the 
EOM to the correct angle on our breadboard system. 

A very slight deviation leads to many other terms in 
the output Stokes vector. The orientation angle calcu- 
lation is very sensitive to these terms. The final flight 
system design will correct this problem by enabling 
rotational alignment of the EOM. 

To get a more accurate measurement of the retarda- 
tion, we placed a precision 7t/8 (0.3927 radians) plate 
in the sample space. This retardation plate has a toler- 
ance of ±0.001 waves or ±0.0063 radians according to 
the manufacturer. The PEM-based system measured 
the retardation at 0.401 radians while the EOM-based 
system measured it at 0.396 radians. 

DISCUSSION AND CONCLUSIONS 

We have developed a stable transverse electro- 
optically modulated birefringence measurement system 
suitable for microgravity fluid physics applications. 

The experimental results from the dual-crystal EOM 
compare favorably with those of the PEM. They also 
agree within the accepted tolerances of the retardation 
plates used. This type of modulation scheme has many 
advantages over existing modulation schemes. Ad- 
vantages include variable modulation frequency from 
0 to 100MHz, low power consumption, low mass, 
small size, and an extremely rugged package. The 
only potential disadvantage is that LiNb0 3 may expe- 
rience photorefractive damage with high-powered la- 
sers. 

Due to its small size, this type of system may be use- 
ful in many types of terrestrial and microgravity appli- 
cations where experiment space is very limited. 
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